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Abstract 

We find analytical solutions to the Cahn-Hilliard equation for the dy- 
namics of an interface in a system with a conserved order parameter 
(Model B). We show that, although steady-state solutions of Model B 
are unphysical in the far-field, they shed light on the local dynamics of 
an interface. Exact solutions are given for a particular class of order- 
parameter potentials, and an expandable integral equation is derived for 
the general case. As well as revealing some generic properties of inter- 
faces moving under condensation or evaporation, the formalism is used to 
investigate two distinct modes of interface propagation in systems with 
a metastable potential well. Given a sufficient transient increase in the 
flux of material onto a condensation nucleus, the normal motion of the 
interface can be disrupted by interfacial unbinding, leading to growth of 
a macroscopic amount of a metastable phase. 

PACS numbers: 64.60.My, 05.07.Ln, 64.60. Qb 

1 Introduction 

The kinetics of phase ordering is a central topic in nonequilibrium statistical physics. 
Much of our understanding is based on theories that describe one or more slowly- 
varying density (or order parameter) variables, governed by a local Langevin equa- 
tion In general, the density variable(s) evolve(s) systematically in response to 
a driving force, which is a derivative of the underlying free energy functional, with 
some mobility (characterized by the Onsager matrix). On top of this are noise 
terms whose magnitude is fixed by requiring that the Boltzmann distribution is a 
stationary state of the dynamics. The nature of the Onsager mobility depends on 
the kind of ordering involved; specifically we must distinguish conserved order pa- 
rameters from nonconserved ones. In the conserved case, the density in some region 
can change only by diffusive transport across its boundary; its time derivative is 
therefore the divergence of a current. This is not the case for nonconserved order 
parameters, which can change locally in direct response to the driving force. 
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The low temperature (noise-free) limit is usually considered appropriate for the 
study of phase ordering kinetics, in which a system is prepared far from equilibrium 
and then allowed to evolve. For example, a uniform high-temperature phase can 
be quenched into a region where it is either locally or globally unstable with re- 
spect to separation into two macroscopic phases. Local instability leads to spinodal 
decomposition 0; if the system is locally stable, phase separation proceeds by a nu- 
cleation and growth mechanism [Q. In either case, the governing equation for phase 
separation of a conserved scalar order parameter is the Cahn-Hilliard equation [Q: 



sion, which treats the order parameter, p, as slowly varying; T is the mobility. Note 
that in principle Cahn-Hilliard theory can accommodate an arbitrary form of /(p), 
which is the free energy density for a homogeneous state. (In particular, it does not 
assume that / is a polynomial in the order parameter p, as would be assumed in 
the time dependent Landau-Ginzburg theory of dynamics close to a critical point 
[|].) Indeed, the approach should be qualitatively applicable even if /(p) consists 
of the lower envelope of several unrelated functions representing phases of different 
symmetry. The free energy near a liquid-solid transition is of this form, for example, 
with p the material density. This assumes only that, whatever other order param- 
eters distinguish the various phases (such as crystallinity), these can relax quickly 
and hence that the rate-limiting process for the phase ordering is transport of p. 
It is conventional, in equation |l|, to treat i^' as a constant (independent of p). We 
do this in what follows, although it might be a dangerous assumption when / is a 
composite function as just described. With this caveat, equation |I] will be relevant 
at long times and large distances if the other order parameters are nonconserved. 

In this paper, we therefore consider the phase ordering problem for relatively 
general forms of /(p), where p is a concentration variable. We assume this is the 
only conserved order parameter, thereby ruling out systems with significant concen- 
tration deviations in more than one species, and also ruling out consideration of heat 
transport. This latter restriction might be severe in (say) metallurgical applications, 
but not for soft condensed matter systems (such as colloidal suspensions) which are 
our main interest. Indeed, for many such systems the latent heats of phase changes 
are entirely negligible 0. 

Our work is motivated by the desire to understand better the role of metastable 
phases in the kinetics of phase separation. That such a role exists has been long 
acknowledged: for example, the "Ostwald rule of stages" |^ asserts that a system 
will progress from an unstable to a stable state, not directly, but by a sequence of 





As explained in Section |], is the square-gradient coefficient in a free energy expan- 
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steps through any intervening metastable states that may be present. In the area 
of metallurgy, there is an extensive folklore on the subject Here we aim at a 
more fundamental understanding, based on a direct analysis of the Cahn-Hilliard 
problem. For the most part, we work in one space dimension. 

Specifically, we shall focus on steady- state solutions of the Cahn-Hilliard equa- 
tion, in which interfaces between phases move with constant velocities. This ap- 
proach appears at first paradoxical, since the diffusive nature of the transport rules 
out true constant- velocity solutions when conserved order parameters are present. 
(Indeed the basic scaling of lengths in diffusive transport is with t^/^ whereas con- 
stant velocities would imply linear scalings.) However, as we discuss later, signifi- 
cant physical insights can be gained by viewing the interfacial motion as a quasi- 
steady process. A broadly comparable analysis, for nonconserved dynamics (where 
true steady-state motion is possible) has been given recently by Bechhoefer et al. 
[^, |IU], |TT|. These authors showed that under sufficient supercooling (or, equiva- 



lently in a ferromagnetic system, sufficient applied field) the interface between two 
stable coexisting phases could become dynamically unstable toward "splitting" . The 
splitting instability results in a macroscopically thick slab of a metastable phase ap- 
pearing between the two stable phases, which can then grow. 



One of the main questions we address here, and in a companion paper |]12 



whether the same scenario is possible for the conserved order parameter case. In 
this paper we show that, although there is no mathematically direct analogue of the 
splitting instability found by Bechhoefer et al., a sufficient transient flux from the 
less dense to the more dense stable phase can indeed cause interfacial splitting. We 
also argue that the split mode will be maintained so long as the supersaturation of 
the less dense phase is sufficiently large. In the companion paper |0 we study the 
long time limit in which the interfaces become sharp on the scale of their separation, 
and give a more detailed discussion of the critical supersaturation required to sustain 
the split mode at long times. That paper also contains a discussion of experimental 
evidence, involving colloid-polymer mixtures [0, which suggests that the onset of 
the split mode might be connected with the observation of arrested crystallization, 
beyond a threshold of supersaturation, in the transition from a colloidal fluid to a 



colloidal crystal. A brief account of these ideas is given in |0 



The rest of this paper is organized as follows. In Section ^ we recall the Cahn- 
Hilliard equation and discuss the conditions under which a one- dimensional treat- 
ment should suffice. In Section ^ we formulate the quasi steady-state form of the 
Cahn-Hilliard equation, paying careful attention to the boundary conditions that 
are required to make the solution physically meaningful. An exact solution is de- 
scribed for a piecewise quadratic potential /(p), with piecewise constant mobility. 
In Section ^ we describe in more detail the properties of the solution, focusing on the 
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case where f{p) shows a metastable minimum of intermediate density. We argue 
that there is no critical velocity above which the interface between stable phases 
ceases to have a steady-state solution (in contrast to the nonconserved case) and 
in Section ^ we show explicitly that such an "unsplit" mode of interface motion is 
linearly stable. In Section however, we show that a split interface, should one 
arise, can also be dynamically stable under appropriate conditions. We discuss qual- 
itatively the nature of the (large) perturbation required to cause splitting. Section ^ 
summarizes our conclusions. Appendix A provides some details of the exact solution 
for the piecewise quadratic case, whereas in Appendix B we derive an exact integral 
representation of the steady-state solution for general potentials, thereby confirming 
and extending some of the earlier results. 



2 The Cahn-Hilliard Equation 

Consider a system characterized by one conserved, scalar order parameter, such as 
local mass density, in a part of the phase diagram where two-phase coexistence is 
the equilibrium state. If the system is far from criticality, and is initially out of 
equilibrium, then its evolution towards equilibrium obeys Model B, described by 
equation |I] (the Cahn-Hilliard equation), which is derived as follows. Let the free 
energy of the system be a functional F[p] of the order parameter p{x). Then the 
chemical potential is defined by the functional derivative, 

6F 
6p 

{n can thereby depend on gradients of p, as well as p itself). Currents in Model B 
are induced by gradients of the chemical potential: 

3 = -TV 11 

where the constant of proportionality, F is the Onsager mobility, which may be a 
function of p. Since the order parameter is conserved, its time-derivative is given by 
a continuity equation, 

dp ^ . 

Let the free energy functional be of the form 

F[p] = jd'x[f{p) + \K{Vpf}. 

The form of the bulk free energy density (or 'order-parameter potential') /, and 
the value of K, are system-dependent. Let the system in question be initially ho- 
mogeneous at a non-equilibrium value of p, between two minima in /; i.e. within 
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a two-phase coexistence region of the phase diagram. For definiteness, let p be ini- 
tially close to the low-density minimum in /. Small fluctuations initially induce 
the evolution of p towards equilibrium. Early stages of the evolution proceed by 
nucleation, if / is convex at the given value of p, or by spinodal decomposition, if 
concave. This stage of the dynamics is not addressed here. By whatever process, 
domain walls (interfaces) soon form. If surface tensions are neglected (legitimate 
when typical interfacial radii of curvature are large, which we assume), subsequent 
motion of a wall is driven by diffusion from the far-fleld. The local profile of the 
wall changes on a shorter time scale than the far-field gradients which determine the 
flux onto the wall, simply because of the difference in length scales: while typical 
distances between interfaces are proportional to t^/^, the characteristic width an 
interface remains of the order of K/ {(Pf/dp^) at all times. Therefore, although 
typical inter-wall distances vary with time, intra-wall dynamics (concerning the lo- 
cal density profile of an interface) soon become approximately steady-state, with a 
quasi-constant input fiux. 

These local interface dynamics are concerned with the movement of a (c? — 1)- 
dimensional wall in a rf-dimensional space, translating normal to itself. Hence, if cur- 
vature and surface-tension effects are ignored, the problem becomes one-dimensional; 
the Cahn-Hilliard equation then reduces to 

at dx \ ^^'dx [dp dxyj' 

In the next Section, we show how to find steady-state solutions of this equation for 
a piecewise-quadratic potential. 

3 Exact Steady-State Solution 

Let us transform to a frame moving at velocity v, in which the position coordinate 
is 

s = x — vt, 

and introduce the notation fn{p) = {d/dp)"'f{p). Demanding that the time deriva- 
tive vanishes in this frame, we find that a steady-state solution p{s) to the Cahn- 
Hilliard equation, travelling at velocity v satisfies 

which is a third-order ordinary differential equation in p'(= dp/ds). If F and /2 
are both independent of p then the right-hand side vanishes and the left-hand side 
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becomes linear and homogeneous in p' . So the problem is (at least) piecewise-soluble 
for a piecewise-quadratic potential /(p), with piecewise-constant mobility F. Such 
a model is actually fairly versatile, and we therefore explore it in detail. 

At discontinuities in our piecewise-constant /2 and/or F, the solutions for the 
separate pieces must respect certain matching conditions. Specifically, the current 
must be continuous so that V ■ j remains finite; the chemical potential must be 
continuous to avoid infinite currents; and continuity of the gradient of p is required 
to avoid infinities in the chemical potential. These three conditions are sufficient to 
fix the constants of integration for the above third-order equation in p'. Subsequent 
integration to find p gives rise to an additional arbitrary constant, which is fixed by 
demanding continuity of p itself. The full solution is given in detail in Appendix 0. 
For each piece of the potential, this solution has the form 

P'(^) = EV^^^ (3) 

where the constants Aj and LOj depend on /2, F and v. 

In addition to the matching conditions (which fix some of the constants arising 
in equation ^ from the division of /(p) and F(p) into pieces) boundary conditions 
are required, to select a specific solution to the differential equation. The correct 
choice and interpretation of these boundary conditions are not trivial; we discuss 
them carefully before proceeding further. The case of interest (depicted in figure [^) is 



pA 




Driving Region Quasi-Steady-State Region 



S 

Figure 1: Schematic diagram of the system. Density p is plotted against distance 
s. On the right of the figure, a high density domain has formed, for which p — ^ p^ 
as s — *• ^oo. On the low density side of the interface, the region for which the 
quasi-steady-state approximation holds meets the 'driving region', which becomes 
depleted as material condenses onto the interface. 

where a region of the high-density phase has formed, and is growing by condensation 
from the supersaturated low-density phase. (Note that, throughout this paper, the 
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high-density phase is depicted on the right of the diagrams, and therefore grows by 
leftward motion of the interface. The opposite convention is adopted in Refs. [|l^ 
and IQ.) The interface is to be modelled in isolation, so the high-density phase 
is semi-infinite. Two boundary conditions arise from this: that the density (and 
hence chemical potential) asymptotes to a constant value, poo as s ^ +00 (say), 
and that the flux asymptotes to zero in this limit. We may now either regard v in 
equation ^ as a given constant and then, from integration, deduce the conditions 
at the other boundary, or v may be seen as an eigenvalue which is set by further 
boundary conditions. In any case, at the second boundary a flux is required (to 
induce motion), and this implies a gradient in the chemical potential. Hence it does 
not make sense to put this boundary at s ^ —00, as fi would be infinite here; instead 
the left-hand boundary must be at some finite position. This raises two possible 
worries: that the thermodynamic limit cannot be taken since the model is of a finite 
system, and that steady-state solutions cannot be found in a finite system. The 
interpretation which resolves these difficulties is as follows. The left-hand boundary 
is at a finite distance from the interface, and moves with the interface. It is not in 
fact the edge of the system, but simply the point at which the behaviour ceases to 
be quasi-steady-state. The part of the system to the left of this boundary, which 
does not solve the steady-state equation, may be referred to as the 'driving region', 
since it is responsible for supplying the quasi-constant fiux and chemical potential 
to the propagating interface, by non-steady-state diffusive depletion of material. In 
summary, a specific solution to the steady-state equation is fixed by the asymptotic 
value of p in the limit s +00 and by the values of the chemical potential and fiux at 
some (rather ill-defined) position to the left of the interface, where the steady-state 
region meets the driving region. 



4 Properties of the Solution 

One important qualitative observation, noticeable on graphs (such as those discussed 
in section |4.1|) of the exact steady-state solution calculated in Appendix ^, is that the 



characteristic width of the interface always decreases as the speed of condensation 
(or equivalently the incident fiux) increases. In Appendix it will be shown, to 
first order in v, that this is a generic result. A second qualitative result, found below 
(Section [4.2|) , is that steady-state solutions exist at all velocities v. This holds even 
when an intermediate metastable well is present in the order parameter potential; 
accordingly (and in contrast to the case of a nonconserved order parameter 0) there 
is no critical velocity above which the interface must split. Both of these results have 
implications for the formation of metastable phases. 
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4.1 Form of the Interfacial Profile 



Before discussing splitting, we show some typical numerical results for a steadily 
moving interface, in a system where a metastable phase is possible. Consider a 
system in which /(p) contains a metastable well at a density between that of the 
growing high-density phase, and the supersaturated low-density phase. A piecewise- 
parabolic form for such a potential is shown in figure |^. ( Met ast ability requires 
that the middle well is above the common tangent to the other two wells.) Let the 





Figure 2: (a) A piecewise-quadratic three- well potential /(p), for which all relevant 
parameters are given. Well number 2 is metastable. (b) The density profile of an 
equilibrium interface for the above potential. Regions are labelled in which the 
density corresponds to wells 1, 2 and 3 of the potential, (c) Density profile of an 
interface moving at velocity —0.027^^ for the same potential, and constant mobility 
r. Notice that region 2 is narrower than in the equilibrium case. 



low, intermediate and high density wells be referred to as 1, 2 and 3 respectively. 
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We may define the order parameter p (chosen as a scaled, relative density) to be ±1 
at the minima corresponding to phases 1 and 3, and take their free energy densities 
/(— l),/(+l) to be equal. (This choice involves no loss of generality, since, as is 
well-known adding a linear term to the free energy density (/ f + ap + b) 
has no effect on the solutions of the Cahn-Hilliard equation.) Steady-state interface 
profiles for this system are given in figure |pD for f = {i.e. the equilibrium wall) and 
(in fig. I^c) for v = — 0.02irr. In each case the mobility F is constant throughout 
the system. Negative-f solutions are of greatest interest, since they describe the 
condensation of material from a supersaturated region onto a growing domain, and 
are therefore central to phase-ordering dynamics. Notice that both solutions shown 
have an inflection at the density of the metastable phase 2. Without forming a 
macroscopic amount of the metastable phase, the interface takes advantage of the 
local minimum in free energy by having extra material at this density. Region 2 is 
noticeably narrower in figure ^ (for condensation) than in figure ^ (the equilibrium 
profile) . 

4.2 Existence of Solutions for an Unsplit Interface for All v 

We are interested in whether such a steady-state interface might split into two parts, 
the 1-2 part of the interface propagating faster than the 2-3 part, analogously to the 
'dynamic splitting instability' ||^ which can arise in the dynamics of a non-conserved 
order parameter in a three-welled potential. In the non-conserved critical 
velocity exists, above which there exists no steady-state solution for the propagation 
of a 1-3 interface; instead a macroscopic amount of the metastable phase 2 must be 
created between a pair of moving (1-2 and 2-3) interfaces. 

In comparing the conserved and non-conserved dynamics however, an important 
distinction should be borne in mind. In the non-conserved case, the velocity of each 
interface is controlled by an external field, which adds a linear term to the potential. 
(Indeed, to obtain a dynamic splitting instability in the non-conserved case, the field 
must cause the potential in the middle well to fall below that of one of the others.) 
In the conserved case, however, linear terms in the potential are irrelevant; instead, 
the velocity is controlled by the boundary conditions. We now present an argument 
showing that the unsplit propagation mode exists for all velocities in this case. (The 
argument is not limited to the case of piecewise quadratic potentials.) 

First, note that equation ^ may in principle be integrated spatially, from right 
to left, for a given v and p^o, to find the value of p at any point. This could fail to 
produce a solution for a three-well potential, only if the resultant profile p{s) fails 
to span all three wells due to the presence of a minimum in the function p{s). This 
would occur whenever the given value of v was above the critical value. However, 
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such turning points in p(s) do not arise in the steady-state solutions. This follows 
from the expression for the chemical potential 

dp dx"^ 

To see why, consider first the equilibrium interface profile, for which /i is a con- 
stant. Clearly, this spans all three wells. Any solution of negative v (describing 
condensation) must have a higher chemical potential than the equilibrium value, 
at any given point on the interface where p < poo, because there is a steady flux 
onto the high-density phase. Hence, for any given value of p, it follows from the 
above expression for p that the curvature of p{s) must be more negative than for the 
equilibrium profile. So no minimum exists in p(s). This argument also holds if poo 
is greater than the equilibrium value, since there is still no minimum to the f = 
solution in this case. The argument may even be extended to the case where poo 
is below the equilibrium value so that the f = solution exhibits a minimum. So 
long as this static solution spans all three wells, negative v solutions with the same 
asymptote must also do so, since their curvature at any given value of p is more 
negative. Hence there is no critical velocity for condensation. 



5 Stability of the Solutions 



Having established the existence of solutions corresponding to unsplit interfacial 
propagation at all velocities f , we now show that these solutions are stable against 
linear perturbations. (The argument below is restricted to piecewise quadratic f{p).) 
Let the field p{x, t) obey the full Cahn-Hilliard equation of motion, and be written 

as 

p{x,t) = po(x,t) + e{x,t) 

where po{x,t) is a solution of the steady-state equation, and e is initially small. 
Differentiating with respect to time gives 

de^ _ 
di " 





PO/ 



where 



f^\po+e = f^lpo + ^/2(po) - K— + 0{e 



For the exact solutions calculated in Appendix 0, /2 and V are both piecewise 
constant. It follows that, on any piece of this solution, a small perturbation e about 
the solution obeys the linearized equation of motion 

— = r/2— - TK— 

dt dx"^ dx"^ 
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In Fourier space, writing 



oo 



e{x, t) 



— CO 



the equation of motion becomes 



de{q, t) 
dt 



{Tf, + TKq')e{q,t). 



In any well in the potential, the coefficient —q^{Tf2 + TKq^) is negative for all q. 
So all Fourier modes of a small perturbation, about any piece of the solution in a 
quadratic well, decay exponentially with time. Also, since the pieces of the solution 
must always obey the matching conditions at cusps in the potential, the solution as a 
whole is stable, when all continuous parts of the potential are convex (as is the case, 
for instance, in the potential of figure ^). Furthermore, in any concave part of the 
potential, solutions are stable with respect to perturbations of shorter wavelength 
than 2'nJKj[—f'2). Hence, so long as the spatial distance over which the interfacial 



proffie spans a concave part of the potential is less than iT^^Kji—j^)^ the solution 
will be stable. It is obviously true that the equilibrium interface profile satisfies 
this linear stability criterion. As stated previously, we have observed solutions, for 
interfaces moving due to condensation, to be narrower than the equilibrium interface 
(and proved it to ffist ffist order in v in Appendix |BI); hence they too are linearly 
stable. 

Although this completes the argument, it is useful to have a conceptual picture 
of the mechanisms giving rise to this stability. Consider once more the proffie in 
figure ^ of a moving interface for a potential with a metastable middle well. A 
perturbation such as that shown by the dashed line in figure ^ may be considered as 
a small change in s at constant p, rather than vice versa. The perturbation shown 
is tending to separate the 1-2 part of the interface from the 2-3 part, thus widening 
the metastable region. Notice that this makes the curvature more negative on the 
part of the interface labelled 'A', and more positive at 'B'. So the chemical potential 
is increased at A and reduced at B, thus enhancing the flux onto the 2-3 part of the 
wall. So the 2-3 part of the wall will catch up with the leading 1-2 part, and steady- 
state motion will be restored. This negative feedback mechanism is responsible for 
the linear stability of an unsplit interface. 

It is worth noting, in addition, that numerical solutions of the Cahn-Hilliard 
equation have been performed, using a three-well potential, and have confirmed 
that 1-3 interfaces may propagate stably, even at large values of v. 
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Figure 3: Density profile of a steadily moving interface (solid line) and the result 
of a perturbation in the direction of the arrows (dashed line) which has increased 
the width of the region occupying the metastable middle well of the potential. The 
result is a more negative curvature in the position labelled A, and more positive 
curvature at B. 



6 Local Dynamics of a Split Interface 

It was shown above that steady motion of a 1-3 interface, spanning the intermediate 
metastable well, is linearly stable. Hence a 1-3 interface, once formed, continues to 
propagate in the absence of large perturbations. Such perturbations may however 
arise, especially in the early stages of interface formation (including the dynamics 
prior to the time at which the Cahn-Hilliard equation becomes a good approxima- 
tion). Therefore, let us consider the situation whereby a large slab of metastable 
phase 2 has formed, by whatever mechanism, so that the 1-2 and 2-3 interfaces are 
separated by a distance large compared with the scale ^Jk/ f2 set by the curvature 
term in the free energy. Such a situation is depicted in figure ^, in which various 
quantities are defined: the interface separation Aa;, the two interface velocities Vi 
and V2, the fluxes into the 1-2 interface, ji, and between the interfaces, j2 (both of 
which are taken to be approximately constant over the spatial regions of interest), 
and four densities pA, Pb, Pc and pn- Let us introduce a further approximation, as 
follows. We assume that the interfaces are moving sufficiently slowly that the den- 
sities {pa, Pb) and (pcPo) on either side of each are approximately the values that 
would arise at coexistence of the two given phases, in the absence of the third. These 
pairs of values of p may be found from the bulk free energy density by the usual 
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Figure 4: Density profile of a system in whicli 1-2 and 2-3 interfaces are separated 
by a distance Ax. The interfaces are travelling with velocities vi and V2 respectively. 
Densities pA, Pb, Pc and po and fluxes ji and j2 are also defined in the figure. 



double-tangent construction, as shown in figure ^ (which also shows the construction 
for the globally stable binodal values {pa, Pp) for 1-3 phase coexistence). In princi- 
ple the double tangent construction is subject to small corrections due to interface 
motion; these are calculated, for completeness, in Appendix |C| but we neglect them 
here. 



6.1 Growth or Collapse? 

The time evolution of such an interface is not strictly a question of steady-state (or 
even quasi-steady-state) dynamics. Accordingly we give only a brief discussion, and 



leave a fuller exploration of this interesting problem to the companion paper [|12 
(see also [Q). The basic issue is whether the slab of metastable phase grows or 
shrinks. 

When curvature is small, as is the case between the interfaces, the Cahn-Hilliard 
equation is well approximated by the diffusion equation, p = DV^p, with the diffu- 
sivity given hj D = T f2{p) which is approximately constant given that p does not 
vary much. {D is exactly constant in a quadratic potential well of fixed F.) Notice in 
figure H that the inequality > Pc is a necessary result of the metastability of well 
2. Hence, since the diffusion equation governs the inter-wall region, j2 is positive. 
Thus flux flows onto the 2-3 interface, contributing positively to V2 and negatively 
to Vi. So, the effect of j2 is to reduce Ax, as would be expected for the dynamics 
of a metastable phase. If phase 2 is to grow, the constant flux ji into the system 
must be sufficiently large to make |fi| > \v2\- Invoking the diffusion approximation. 
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Figure 5: Figure demonstrating the double-tangent construction to find the stable 
binodal densities (po,, pp) and the metastable binodal densities (pA, Pb) and (pc, Pd) 
in a three-well potential. This potential can be converted to the form of fig.0 by 
adding or subtracting a linear term (/ f + ap + b), and then both shifting the 
origin of p, and rescaling it. 



and the linearity of the function p{x) in phase 2, the condition for growth of the 
metastable phase (rather than recombination of the 1-3 interface) becomes 

,.>fl(££_Mfi + ^iL^). (4) 
Ax V pD- Pc) 



When this condition is satisfied, a 'split' mode of interfacial propagation can 
arise, which is fundamentally different from the propagation of an unsplit 1-3 inter- 
face in a number of respects. The most important distinction is that now, if ji is held 
constant, the width of the metastable region. Ax grows without limit. In contrast, 
in a stable 1-3 interface. Ax remains bounded. Indeed, in the equilibrium {v = 0) 
interface. Ax is of the order of the characteristic interfacial width ~ ^K/ f2, and 
at higher speeds. Ax becomes smaller. Accordingly there is an upper bound on Ax 
close to the equilibrium value (although not equal to it, since unsplit propagation is 
resumed after a small, positive perturbation in Ax). 

Another qualitative difference between the split and unsplit modes of propagation 
is their response to a perturbation. It was demonstrated in Section ^ that increasing 
the width Ax of the metastable region led to an increase in the flux though it, 
resulting in a negative feedback mechanism. On the other hand, if the interface is 
split, and therefore non-monotonic and containing a well-developed region in which 
diffusive motion is dominant over curvature-induced motion, increasing Ax reduces 
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the gradient in region 2. This reduces j2, and causes the 2-3 interface to lag still 
further behind the 1-2 wall. So the corresponding feedback in split interface motion 
is positive. It follows that, at constant ji, there is a barrier (in configuration space) 
to the formation of an unsplit interface, but once this barrier is crossed, such an 
interface will remain split indefinitely. 



6.2 Selection of Split or Unsplit Mode 

It has been shown that (at least for piecewise-quadratic potentials) the propagation 
of an unsplit interface is locally stable, and that, given sufficient input flux, the split 
interface mode is also "stable" (in the sense of remaining split indefinitely). The 
question arises of which mode of evolution will be selected in a given system, and 
how it might be possible to change from one to the other. Clearly in a real system, 
the fiux input to an interface is not constant. Normally, in late-stage evolution, it 
is a diminishing function of time. One might conclude from this that the criterion 
for growth of the split mode (equation H) must at some point be violated. However, 
this criterion becomes easier to satisfy as Ax increases. The ultimate fate of a split 
interface in fact depends on supersaturation: this is described in the companion 
paper |1^ . 



The problem of how to 'unbind' a 1-3 interface is somewhat clearer: a large 
transient increase in input flux is required, to overcome the negative feedback mech- 
anism described in Section ^. There is presumably some critical value of Ax at 
which the feedback switches from negative to positive and the interface splits. The 
transient increase in flux must be sufficient to separate the leading (1-2) part of 
the interface from the trailing part by this critical amount, before the interface can 
deliver a restorative increase in flux to region 3 by curvature-induced motion. If 
the transient increase of input flux is insufficient, and the 1-3 interface adjusts to 
the new higher speed, the criterion for unbinding it becomes more stringent (since 
higher-speed interfaces are narrower, and hence both further from the critical value 
of Ax, and 'stiffer' in terms of the negative feedback mechanism). Since fluxes tend 
to decrease with time during late-stage phase-ordering, a transient increase in flux, 
sufficient to unbind a 1-3 interface, is most hkely to occur during the early stage 
dynamics (nucleation or spinodal decomposition). These dynamics are not quasi- 
steady-state, and we do not discuss them further here. But it is interesting that, 
whenever metastable phases are possible, the details of these early stage dynamics 
can determine the gross features (split vs. unsplit mode) of phase separation at 
much later times. 

Finally, in the context of mode selection, a useful distinction can be drawn be- 
tween two types of interfacial binding/unbinding. The type described above can be 
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called "curvature unbinding" - the process whereby Ax becomes large compared to 



K//2. In the companion paper |T2[ (see also |jTj]) we study "diffusive unbinding" 
which is hnked to the evolution of ji (treated above as an externally imposed pa- 
rameter). In Ref. |]14|, we also report numerical results which show that curvature 



unbinding certainly does occur, within the Cahn-Hilliard equation, at least for some 
parameter values and some initial conditions. These include cases (such as an initial 
step function wall between phases 1 and 3) which, though strongly perturbed from 
the equilibrium profile, are definitely not unbound to begin with. 

7 Conclusion 

It is often assumed that steady-state solutions of the Cahn-Hilliard equation (Model 
B), for the phase ordering dynamics of a conserved order parameter, are unphysi- 
cal and therefore uninteresting. In this article we have shown that such solutions 
shed light on the local dynamics of interfaces, so long as care is taken to inter- 
pret the boundary conditions correctly. Exact solutions were found for the case of 
any piecewise-quadratic order parameter potential and piecewise-constant mobility; 
these are locally stable. For the most general case of arbitrary bulk free energy den- 
sity and mobility as functions of mass density p, a systematic expansion scheme was 
derived (Appendix to solve Model B for a moving interface. Using this scheme, 
it can be shown (to at least first order in v) that an interface contracts when moving 
under condensation, and expands during evaporation. If a metastable phase exists, 
whose density is intermediate to the two stable phases, a moving interface between 
the stable phases is still locally stable with respect to small perturbations at all 
speeds, but may be 'split' by a large, transient increase in the flux of condensing 
material, and thereafter exhibit a qualitatively different mode of propagation. Such 
a transient disturbance is most likely to arise in the early stages of nucleation and 
growth. The split mode of interface motion results in the formation of a macroscop- 
ically large amount of the metastable phase and relies on a sufficiently large flux 
of condensing material being maintained. A companion paper (see also [p!^). 



discusses further the implications of these findings for the growth rates of competing 
stable and metastable domains during the phase-ordering process. 
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A Exact Solution of the Steady-State Equation 
for Piecewise-Quadratic Potential 

Equation H is now solved to find p(s) for any piecewise-quadratic potential /(p) and 
piecewise-constant mobility r(p). Between discontinuities in /i, /2 or F, p(s) solves 
the equation 

//// // I / n 

p — ap ~ op = 
where a = f2/K and b = v/ {KT). The solution is 

for some constants Aj, fixed by the boundary conditions and matching conditions. 
The constants Uj are the three roots of the cubic equation 

uj^ — aujj — 6 = 

which is solved, for 276^ < 4a^, by 

^ fa . / a — 271 
^1 = 2^-sm(^^ 

fa ( o 

fa {a + 2-K 
= 2^-sin(^^ 

where 

cr = — arcsm 



4a3 ■ 

If 276^ > 4a^, the solutions may be written in terms of the real quantities 

a = {\h + yb^-±a^Y/^ and f3 = {^b - i^fe^ - ^^3)i/3 

LUi = a + (3 

Us = UJ2 

Ml \ 

All that remains to be found is the vector of coefficients A = A2 \ for each 
piece. This is fixed by assigning values to p', p" and p'" at the right-hand boundary 
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of each section, which may be re-defined as the origin of s by multiplying A by an 
appropriate exponential factor. From the equation for p'{s), it follows that 




/111 

A = UJ\ UJ2 

\ ujf Ji 

The vector of derivatives of p is determined, from the solution in the neighbouring 
piece, by the matching conditions given in Section ^ where it was stated that 
the gradient, chemical potential and fiux must all be continuous. That is, the 
quantities p', [Kfl' — /i) and {K(^" — f2p')^ are continuous. These quantities can 
be calculated from the neighbouring solution, given the value of s at which it meets 
the discontinuity (say s = Si). If the discontinuity (in the mobility or potential) 
is at a value p = pi, then si is given by inverting the equation p(si) = pi which, 
unfortunately cannot in general be done analytically. (However, the inversion is 
always possible for the reasons explained in Section |4.2i ) Hence one numerical step is 
required in the solution. (Of course the special case of a double-well potential, with 
just one cusp discontinuity, is completely soluble analytically, since the arbitrary 
origin of s can be put at the cusp of /(p).) 

Finally, the boundary condition that p — > poo as s — +00 results in the solution 

/ -(Poo - Pi) 

for the highest density section: ^ = | and hence 

V 

P'is) = -(Poo - Pi) ^^3 e'^''* for s > 0, and p(0) = pi. 



B Integral Solution of Cahn-Hilliard 

In Section ^, the assertion was made that interfaces generically contract when mov- 
ing under condensation. This result and its converse (that interfaces expand during 
evaporation) will be derived now, using a systematic expansion of steady-state solu- 
tions of the Cahn-Hilliard equation in powers of the interface velocity v. As previ- 
ously, we consider the one-dimensional Cahn-Hilliard equation in a frame moving at 
velocity v, in which distance is measured by the coordinate s; the order parameter 
p asymptotes to a finite constant poo as s — +00. However, in contrast to the 
analysis in Section ^ and Appendix 0, no particular form will be assumed for the 
order parameter potential /(p), and likewise the mobility r(p) will be an arbitrary 
function. 
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B.l Inversion of Variables 

It is convenient to invert the equation so that p becomes the independent variable, 
and the equation is solved for s. This means that the unspecified potential /(p) 
and mobility r(p) are now functions of the independent variable. Henceforth let the 
curvature constant K be set to unity without loss of generality. (This is equivalent 
to measuring time in units of K and length in units of \/K.) After inversion, the 
full, one-dimensional Cahn-Hilliard equation in a moving frame may be written 

ds d 
dt dp 

Clearly (as is well-known) adding a linear term to the potential f{p) has no effect. 

B.2 Integration of the Steady- State Equation 

Let us now set the time derivative to zero in this moving frame, and integrate once 
with respect to p, to obtain the steady-state equation 

where the arbitrary linear term — (Ap + c) has been explicitly added to /(p). Note 
that the constant A is equal to the chemical potential in the asymptotically fiat 
region: A = p(oo) (remember p is non-uniform in a moving interface). It can be 
confirmed that — f Poo is the correct constant of integration, given that the left-hand 
side of the above equation is simply the fiux j. 

Let us define h{p) = ds/dp, so that 1/h = Vp. Then integrating with respect to 
p, we find 

A(i/,-2 _f + xp + c)=v [ ^P-P^^^ dp + const. 
dp J T 

It is easy to confirm that the left-hand side of this equation is the chemical potential, 
measured with respect to the value at poo- Now consider the factors in the integrand. 
As p ^ Poo, {p — Poo) tends to zero linearly, while it is expected (and confirmed in 
the special case of appendix]^) that h{p) —>■ oo logarithmically. Hence the integrand 
vanishes as p ^ poo- So the constant in the above equation is zero, and we have 

-{,h -f + Xp + c) = -vJ^ -^^h{p)dp. 

Let us henceforth absorb the terms Ap + c into the definition of /(p). (That is, 
the arbitrary linear part of /(p) is defined so that the potential has zero value and 





-/ 
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gradient at poo-) Integrating by parts once gives the final result 

""^^TiTwr /wi — W) — 

Notice that the Cahn-Hilliard equation, in its differential form, is fourth order, 
but that this integral representation of the steady-state equation contains only one 
integration. 

If h{p) is written as a power-series in f, then the binomial in equation ^may be 
expanded, and the formula iterated to produce a systematic series approximation 
for h, to arbitrarily high powers of the velocity. The power-series expansion of the 
steady-state solution s(p) is then obtained by integration as: 



s(p) = j Hp) dp- 

Clearly the origin of s (and hence the constant of integration) is arbitrary. 

B.3 First-Order Correction to Interfacial Width 

This method is now applied to find the first-order correction to the width As of 
a moving interface, defined as the distance between two points on the interface at 
which the densities have certain fixed values pi and p2 (which could be the densities 
of the two maxima in a three- well potential, for example). It transpires that 



As = Aso + CV + O(w^) 
P2 dp 



where Asq — , ; 

^V{P) 

and c = ij^^dpfipy-' dp'{p^-p'){p' -p)T{pr'f{p')-'. 

This last equation expresses the rate of change of width of the interface with velocity, 
c in terms only of the two functions which characterize the physics of the system; the 
bulk free energy density f{p) and the mobility r(p). All factors in the integrands 
of this expression are positive over the ranges of integration, so c is positive. Neg- 
ative values of v correspond to growth of the asymptotically fiat, dense region by 
condensation, and positive values correspond to evaporation. So, without assigning 
any special properties to the functions / and F, it has been shown that the interface 
contracts (As < Asq) during condensation and expands during evaporation. 

As mentioned at the end of Section ^, poo may vary from its equilibrium value 
{i.e. the density at which phase 3 coexists in equilibrium with phase 1 for a three- 
well potential). In the present section, a formalism has been developed by which 
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the solution for a moving interface is expanded about the v = solution. Does 
this allow for variation in poo? The answer is 'yes' because there is in fact a whole 
family of f = solutions with different values of poo, of which the equilibrium 
solution is just one member. The equilibrium solution is the special member of this 
family for which p asymptotes to a finite constant as s — oo, rather than growing 
exponentially (either positive or negative) and thus remaining curved and having 
unphysical boundary conditions at s ^ —oo. But, for our purposes, the whole 
family may be used since, as discussed in Section |^, the left-hand boundary is not 
put at negative infinity. 

C Correction to the Double- Tangent Construc- 
tion for a Moving Interface 

The double-tangent construction illustrated in figure |^ gives the equilibrium densities 
of two coexisting phases. Recall its elementary derivation as follows |^ . In a uniform 
part of the system of volume V, containing N particles, the local density is p = N/V, 
and the free energy is F = fV. From these two simple relations, and the definitions 
of chemical potential p = {dF/dN)y and pressure P = —{dF/dV)^, it follows that 
f = fip — P, which is the equation of a straight line on a plot of / versus p, with 
gradient p = /'(p) and intercept —P. Given that two coexisting phases have equal 
chemical potentials and pressures, it follows that the straight lines tangent to /(p) 
at the respective coexistent densities have equal gradients and intercepts. Hence 
they are the same line. 

Consider now the non-equilibrium case of an interface in uniform motion at 
velocity v. By continuity, the flux at any point on the interface is j = — (poo — p)v- 
In model B, Vp = — j/r. Hence, integrating across the interface, the difference 
in chemical potential between the dense, asymptotically uniform phase at positive 
infinity, and any given point A is 



from which the equilibrium result follows for v = 0. 

Now let us consider the local quantity P = pp — f (which reduces to the usual 
definition of pressure in a homogeneous system). With this definition. 




dx 



VP = pVp - KV'^pVp. 



Using again the above expression for Vp, and integrating gives 
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If the curvature part of the chemical potential is small at point A, this expression 
for the pressure difference becomes 



J 



A r 



So the equality of chemical potentials and pressures across a static interface, which 
gives rise to the double-tangent construction, has the following corrections, to first 
order in v, for a moving interface: 

(Poo - p) 



An 
AP 



dx 



°° (Poo - P)P 



dx. 



lA V 

Notice that both corrections depend on the range of integration, i.e. on the position 
of the point A. This is no surprise, since there must be gradients in the pressure 
and chemical potential outside the interface in order induce motion. Notice also 
that both coefficients of v are positive, so Aji and AP are of the same sign as v. 
This is negative for condensation and positive for evaporation. So, as expected, 
during condensation the pressure difference and chemical potential difference across 
an interface are smaller than at equilibrium. 
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